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1. INTRODUCTION 

This paper presents a general methodology for estimating a stratum's 
(observable) crop acreage proportion for a specified season of a target year 
from the crop's estimated acreage proportion for sample segments from within 
the stratum. Sample segment data for several years and seasons are generally 
used in conjunction with those for the target year and season. The proposed 
methodology is an application of estimation from a mixed Analysis of Variance 
(AOV) model . 

The specific model proposed is described and developed in section 2 of this 
document. A general discussion of estimation using the model is also presented 
in this section. In section 3, the use of the model in crop acreage proportion 
estimation is described. The general applicability of the development is 
illustrated by three examples. A documentation of the Statistical Analysis 
System (SAS) implementation of the methodology and sample runs for the examples 
of section 3 are presented in section 4. 

In the following discussion, matrix notation is used extensively. In particu- 
lar, capital letters (with or without subscripts) refer to matrices. Under- 
scored small and greek letters (with or without subscripts) refer to column 
vectors. Small and greek letters that are used with subscripts and without 
underscores refer to elements cf vectors and matrices. Constants are denoted as 
small or greek letters with neither subscripts nor underscores. The transpose 
of a matrix, say A, is denoted A'. The n-by-n identity matrix is indicated by 
Ip. If A is a square matrix. A" denotes any matrix such that AA"A = A. For A, 
n-by-n nonsingular, A“^ denotes a matrix such that AA"^ = A"^A = Ip. The 
symbols ^ and 2p represent column vectors of length n composed entirely of 
zero's and one's, respectively. When the dimensions of a matrix or vector are 
not specifically noted, it is assumed that they are conformable with the 
designated operations. Mathematical expectation is denoted E[ ] for both matrix 
and scalar operands. Finally, if f( ) is a function and is a vector the 
elements of which are denoted p^ , then £(£) denotes a vector such that f^ * 
l'(Pi)- 
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THE BASIC MODEL AND ITS USE IN ESTIMATION 


2 « 


The proposed methodology is an application of the mixed AOV model; 

^ - Xs Zr + e (2.1) 


where 

^ « an n-by-1 response vector 

■ a t-by-1 vector of fixed-effects coefficients 
X = an n-by-t matrix representing the fixed-effects design 
jr* an s-t -1 vector of random-effect coefficients 
Z • an n-by-s matrix representing the random-effect design 
^ » an n-by-1 vector of unexplained errors 

The usual assumptions regarding this model are as follows: 

a. r_ is a random vector with mean 0^ and variance covariance (V-C) matrix V^. 

b. ^ is a random vector with mean ^ and V-C matrix Vg* 

c. and ^ are independent. 

The particular application considered in this document has the property that 
is a transformation of an observable vector that is, = f(£) for some 
function f(»)* Hence, in the following paragraphs, £is referred to as the 
response vector and is referred to as the transformed- response vector. The 
elements of the response vector £ are assigned to have the following properties 
for each i(i * 1, 2, •••, n) : 

d. 0 < p^ < 1. 

e. E[p^] - 

f. Var[p^] « e[(p^ - « n^(l - 71.). 
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The following additional properties are assumed also. 

g. Vg is a diagonal matrix with the ii^i^ diagonal entry proportional to 
0 ^ Var[f(p^)] for some > 0. 

h. Vp is a diagonal matrix with the i^^ diagonal entry equal to for some 

Y > 0. 

i. X and Z are both of full column rank [although (X:Z) need not be of full j 

1 

coluiTw rank]. 3 

j. n > rank (X: Z) + 1. 

The choice of a transformation function f(») will also be restricted to one of 
the three forms: 

k. 'i dentity' , f(p^ ) = Pi . 

l. 'log', f(Pi) = tn(pi). 

m. 'logit', f(Pi) = .5 X tn[pi/(l - Pi)]. 

2.1 DEVELOPMENT OF THE ESTIMATION EQUATIONS 

Under the conditions imposed above, the mean vector and V-C matrix, 
respectively, for ^ are 

E[x] = X0 = ji (2.2) 

and 

Var[^] = E[(^ - ji)(y. - i)'] - Vg + ZV^Z' (2.3) 

-12 2 

By denoting = W and V^ = IgOgY* equation (2.3) becomes 

VarCy.] = (W^ + yZZ' ) o^ * (2.4) 

Hence, the least squares estimator of jB (given W and y) is 

S = (X'VX)‘^XVy^ (2.5) 
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r 


and its V-C matrix is 


Var i 


(X'VX)'^Og 


( 2 . 6 ) 


If C is a q-by-t matr'x of estimable contrasts, then an unbiased estimate of 


iLc - Ci 


(2.7) 


and 


V»''[y - C(X'V)()‘'c'o'- • 


( 2 . 8 ) 


Since the quanti^-v^s W, y, and are usually unknown, equations (2.7) and 
(2.M) can or’y c: for estimation fran tfe mixed model (2.1) after sub- 


A X 


stitut'*',- ,)t their respective estimates, W, f, and o^. (For a discussion of 
the e^^ioct 01 thi? substitution on the asymptotic properties of the least 
squares estimator of shown in equation (2.7), see reference 1.) 

2.1.1 ESTIMATION OF W, o^, AND y 

Note that W = V“^a^ where the i^*’ diagonal entry of V is 
e e * 


Mi 


hence, W ■ [w.jj] where 


Mj 


Var[f(p^ )] • Of(M^®e 


lly/o2(w^) if i - j 
fO if i # j 


(2.9) 


( 2 . 10 ) 


Since E[fJ ■ 3 initial estimate of W can be obtained by replacing each 
in equation (2.10) by its respective estimate p^. 
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That Is, take w > 


P.] 


Mhere 


'U 


i/4(p,) 


if 1 • J 


0 if 1 # j 


( 2 . 11 ) 


If equation (2.1) is premul tipi led by where M ■ the result is 


- i /2 - 1/2 - 1/2 - l /2 - 1 / 2 , . / J \ 

W ^-W Xa + W 2r^-W e-W (X:Z)I”l+e^ 


(i. 12) 


A11 the assumptions of model (2.1) hold for equation (2.12) so that now 

Var[e ] - - I (2.13) 

-w e n c 

Therefore, an initial estimate of is the residual mean square after fitting 
j8 and £ in equation (2.12) 

\ ■ (X:Z)C(X:Z)'W(X:Z)]‘(X:Z)'wj£/[n - rank(X:Z)] (2.14) 


which is unbiased for o when W « W. 

e 


The estimates of W and can be refined iteratively as follows. First, use 
the estimates of j and £ from equation (2.12) to determine an estimated 

A 

transformed-response vector ^ for the model, that is 


£- (X:Z)[(X:Z)'W(X:Z)3‘(X:Z)'W;t (2.15) 


Then define £ " 



where 



(2.16) 


A 

After replacing each p^ in equation (2.11) with its corresponding p^ from 
equation (2.16), the entire procedure can be repeated to produce refined 
estimates of W and o‘. The iteration process can be continued until W (and 
hence o^) "stabilizes." In the following paragraphs, W and o? denote the 
final stable estimates of W and o^, respectively. 
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A A O 

After w and o are available, y can be estimated via the variance component 
c 

analysis procedure known as Henderson's Method 3 (ref. 2). The resulting 
estimator for y Is 


- ( 


y.'W{(X:Z)i:(X:Z)'W(X:2)]’(X;2)’ - X(X' WX)"^X' 


'■/"I) 


rank (X:Z) + s 


tracejz'w[l^ - X(X' WX)‘^X'w]zj 


(2.17) 


which Is unbiased for y when W « W. 


AO ^2 

The estimator in equation (2.14) has the property that > 0 as required. 
However, the estimator y In equation (2.17) may be either positive or nega- 

A 

tive. In the SAS In^lementatlon of this methodology, any y < 0 results In the 
termination of processing, since theoretically y > 0. 


2.1.2 INVERSE ESTIMATION 

Recall that represents a transformation of an observable response vector, _g.. 
Hence the primary Interest in estimation may not be 

(2.18) 

but 

a-r^(i:^) (2.19) 

A straightforward estimator of ^whlch employs equation (2.7) is 

i ■ r‘(ic) 

A 

since f(*) may not be linear, the estimator £ In equation (2.20) may be biased 
for and Its exact matrix of mean squares (ftlS) may be very difficult to 

A 

determine. However, reasonable estimates of the bias and ^MS for £ can be 
obtained as follows. [This development parallels that shown by SI el ken and 
Dahn (ref. 3)]. 
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Consider the approximation of each a 


series expanded about so that 




as a third order Taylor 




(^'i) 


( 2 . 21 ) 


E y 1 ■ y. , taking the expectation of both sides of equation (2.21) yields 

L ‘^ij ^1 

e[.*,] . . (1/2) V.r ^ * .5 . 


SO that 


bias 




a. • .5 K V f 
^ «=11 


(-2)/: 




(2.23) 


The matrix representation Is 

V 


bias 


[i] 


■11 


0 ... 0 

V 

^22 




(2.24) 


qq 

A 

To determine an estimate of MMS for drop the third term froin equation (2.21) 
so that 


a^ ■ 


) 
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(2.26) 

2.1.3 RESPONSE TRANSFORMATION EQUATIONS 

All equations presented thus far are valid regardless of the transformation 
employed to get from £ to provided all the quantities necessary for the 
estimation are defined. For the crop acreage proportion estimation problem, 
three transformations are potentially useful; the "Identity," the "log," and 
the "logit." The formulas associated with these transformations and the 
peculiarities associated with their use are presented In sections 2. 1.3.1, 

2. 1.3. 2, and 2. 1.3. 3. For a more complete treatment of the development of 
these transformations, see reference 3. 

2. 1.3.1 The Identity Transformation 
The Identity transformation Is: 

f(p^)-P^ (2.27) 

From the previous development. It follows that the W matrix associated wi h 
equation (2.27) has diagonal elements 

- 1/»^(1 - w^) (2.28) 

In order to avoid numerical problems in the computation of w, the SAS 1mp1emen> 
tatlon replaces any p^ < cz by ez and any Pi > 1 - c2 by 1 > tz '«'hen evaluating 
equation (2.U). The parameter cz Is specified by the user. (The substitution 
Is also performed, If necessary, during each reweighting Iteration.) 




Since equation (2.27) is a linear function, it also follows that 


A A 


1= ^C 


bias 


[i]- 

MMS [a] = Var [a] = 


0 

"q 


(2.29) 


2. 1.3. 2 The Log Transformation 
For the log transformation 

f(p.) = An(p.) 

This implies the associated W matrix has diagonal entries 

w-i = ir./(l - ir.) (2.31) 


(2.30) 


Also, since equation (2.30) is not linear, its application to equations (2.24) 
and (2.26) yields 




bias 


[i] = (1/2) 


'11 


0 

V 


'22 


• • 


qq 




(2.32) 


MMS 


[i] 


e'^c, 


e^c. 


e^c. 


e'c 


1 ••• 


• * 


y 

e-'c. 
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In the SAS implementation, several steps are taken to avoid numerical 
problems. First, the calculation of W is handled as it is in the case of the 
identity transformation. Second, since iin(p^) is not defined for p^ - 0, a 
"working log" transformation is used which defines the elements of as 

/tn(p^) if p^ > ej 

- ) (2.33) 

Hr) " ^ 0 . P, < 

The effect of this working log transformation is illustrated in figure 2-1. 
(The parameter ej. inust oe specified by the user.) Third, to counteract the 
effect of using equation (2.33) rather than equation (2.27), the elements of 

A 

a are taken to be 


e^Ci if y^. > itn(e^) 




in 




A 

^1 ^c. 


(2.34) 


which is the inverse of equation (2.23). 


2. 1.3.3 The Logit Transformation 
For the logit transformation 


(A) 


f(p^) - (l/2)m 


and the associated W matrix has diagonal entries 

^11 ■ n^(l - it^) 


(2.35) 


(2.36) 
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A1 so. 


bias 





I 

o 

• 

• 

• 

o 

u 

> 

I 

2e^'c,(l - 

fa1- 

• 

CM 

CM 

U 

> 

o 

LiJ - 

I 

o • • • 

• 

• 

• 

• 

< • • 

o 

if 

(l * 


MSS 



where 








1 + e 



(2.37) 


As before, the SAS implementation avoids difficulties in computing W by 
replacing values of p^ less than C 2 with C 2 (greater than 1 - eg with 1 - eg)* 
Also, since equation (2.35) is not defined for e{0, 1}, a working logit is 
used which defines the elements of ^ as 


^i 


P :rr + -m r if 0 < p. < e, 

r -il ^(^ ^ 

je^ + (1 - ep exp|(l - e^) ^jj 

< 1 - 


(2.38) 


-e. 


(P^ - 1) 


1 F :i=jT" 2e,d - e) 

je^ + (1 - cj) exp (1 - cj) ^ j 


if 1 - ej < p. < 1 




The effect of this working transformation is illustrated in figure 2-2. In 

A 

order to counteract the effect of using equation (2.38), the elements of £ are 
taken to be 



(2.39) 
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1 



3. CROP ACREAGE PROPORTION ESTIMATION USING THE MIXED AOV MODEL 


The development of the proposed methodology has been of a very general nature. 
Consequently, it has very broad applicability; and, as will be seen In the 
following paragraphs, the interpretation of the estimation results can vary 
widely depending on the data used and the design matrices (X,Z) specified. 

The remainder of this section will be devoted to presenting the analyses of 
specific models and data. All of the examples will be derived using subsets 
of the data shown in table 3-1. 


Example 1 : Suppose an estimate of the stratum's at-harvest crop acreage 
proportion is desired for each year having representative segment data in 
table 1. The model to be used is 


ijk 




a. 

1 


6 , 

J 


+ r, + e. 


i jk 


where 


(3.1) 


y^jl^ » the transformed crop acreage proportion estimate for segment k at 
growth stage j in year i 

0 | ■ the stratum's transformed crop acreage proportion in year i 

* an adjustment for the tendency of the crop acreage proportion estimate 

J 

for growth stage j to be different from the stratum's at-harvest crop 
acreage proportion 

r|^ « an adjustment for the tendency of the crop acreage proportion estimate 

for segment k to be different from the stratum's crop acreage 
proportion 

e^jl^ « all other unexplained effects 

The year and growth stage effects are fixed; the segment effect is random. 
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TABLE 3-1.- UATA FOR CROP ACREAGE PROPORTION ESTIMATION EXAMPLES 


Observation 
number ' 


Segment crop 
acreage proportion 
estimate 


0.279 

.154 

.149 


Ic'entification of — 


Growth 

Crop year Segment Substratum 




As estimated from Landsat multi spectral data 
h = Midseason; 2 * at-harvest 
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The matrix representation of equation (3.1) using the data shown 1n table Z-1 
is 


in(.279) 

in(.154) 

tn(.149) 

ln(.074) 

nn(.073) 

in(.229) 

tn(.212) 

4n(.275) 

in(.152) 

ftn(.073) 

ln(.069) 

_«n(.210) 


10 0 0 
0 10 0 
0 0 10 

10 0 0 
0 10 0 
0 10 0 
0 0 10 
10 0 1 
0 10 1 
10 0 1 
0 10 1 
0 0 11 


1 0 

0 1 
0 1 



0 0 
0 0 



0 0 
0 0 
1 0 


0 1 
0 0 


0 0 
0 0 


0 

0 

0 

1 

1 

0 

0 

0 

0 

1 

1 

0 


0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

1 



®121 


«222 


®322 


®123 


"a" 


6223 


^2 


6224 


>“3 


®324 




6111 


®212 


®113 


®213 

J 

_®314_ 


(3.2) 


Let C be defined as 


C 



0 0 0 

1 0 0 

0 1 0 

-10 0 
0 1 1 


(3.3) 


A 

Then Is a 5-by-l vector, and Its elements are Interpreted as follows: 

A 

a. ■ an estimate of the stratum's at-harvest crop acreage proportion for 
year 1. 

A 

b. a^ ■ an estimate of the stratim's at-harvest crop acreage proportion for 
year 2. 

c. a^ ■ an estimate of the stratun's at-harvest crop acreage proportion for 
year 3. 
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A 

d. > an estimate of the ratio of the stratun's at-harvest crop acreage 
proportions for years 1 and 2. 

e. ag ■ an estimate of the stratum's observable crop acreage proportion 
during growth stage 1 of year 3. 

Note also that 

A A 

a. a^ - d 2 is an estimate of the change in the stratum's at*harvest crop 
acreage proportion from year 2 to year 3. 

A A 

b. a^ - is an estimate of the change in the stratum's observable crop 
acreage proportion from growth stage 1 to growth stage 2. 

Since the SAS implementation provides the entire matrix of mean squared errors 

A 

for a, measures of dispersion can be calculated for all the estimators listed. 
Example 2 ; Suppose the model is now 

* 0, * ‘ '•k ♦ 'ljkk 

wnere all terms are as defined for equation (3,1), with 

a. Pjj « an adjustment for the tendency of crop acreage proportion estimates 

in substratum i to be different from the stratum's crop acreage 
proportion. 

b. ■ an adjustment due to the interaction of fij and p^. 
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Since X must be a fuH rank, a reasonable matrix representation of equation 
(3.4) using the data of table 3-1 is 

\n(.279) 

£n( .154) 

£n(.149) 
in(.074) 
tn(.073) 
in{.229) 
in(.212) 

£n(.275) 
in(.i52) 

£n(.073) 

Jin(.069) 
in(.210) 


~1 

0 

0 

0 

1 

0 " 


~1 

0 

0 

0 ’ 


0 

1 

0 

0 

1 

0 


0 

1 

0 

0 


0 

0 

1 

0 

1 

0 


0 

1 

0 

0 


1 

0 

0 

0 

-1 

0 

”“1 


0 

0 

1 

0 


0 

1 

0 

0 

-1 

0 

»2 


0 

0 

1 

0 


0 

1 

0 

0 

-1 

0 

“3 

•f 

0 

0 

0 

1 

f'Z 

0 

0 

1 

0 

-1 

0 

h 


0 

0 

0 

1 


1 

0 

0 

1 

1 

1 

n 


1 

0 

0 

0 

U4_ 

0 

1 

0 

1 

1 

1 

\jon^ 


0 

1 

0 

0 


1 

0 

0 

1 

-1 

-1 


0 

0 

1 

0 


0 

1 

0 

1 

-1 

-1 


0 

0 

1 

0 


0 

0 

1 

1 

-1 

-1 


0 

0 

0 

1 



®1211 
®2221 
®322i 
®1231 
®2232 
®2242 
®3242 
eilll 
®2121 
®1132 
®2132 
®3142 
(3.5) 


Let C be defined as 


(3.6) 


1 0 0 0 0 0 
0 10 0 10 
0 0 10-10 

then a is interpreted as follows: 

• an estimate of the stratun's at-harvest crop acreage proportion for 
year 1. 


a. a 


1 


b. I 2 " if' estimate of the at-harvest crop acreage proportion for substratum 
1 in year 2. 


c. a - an estimate of the at-harvest crop acreage proportion for substratun 
3 


2 in year 3. 
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Example 3 ; Suppose an estimate of the stratum's at-harvest crop acreage 
proportion for year 2 which uses only data from growth stage 2 In year 2 Is 
desired. 


An appropriate mixed AOV model Is 

J'k ■ * ''k * *k 

The corresponding matrix representation of the data Is 


(l/2)tn 



1 


1 

0 

0 


'*1 



(l/2)tn 

f.073\ 

- 

1 

u + 

0 

1 

0 




®2 

(l/2)tn 

(.229 J 
1.771 /_ 


1 


0 

0 

1 


/3_ 


_®3_ 


(3.7) 


(3.8) 


This representation, however, does not satisfy the requirement that 
n > rank(X:Z) + 1. The requirements for estimation can be met If the 
fixed- effect model 

yk " u + (3.9) 


Is considered instead. 

The removal of the random effect from the model causes no computational diffi- 
culty— the estimation of y need only be bypassed. In addition, if the true 
(unknown) value of y is small relative to (also unknown), the models (3.7) 
and (3.9) are roughly equivalent in their effect. In the SAS implementation 
of the methodology, specification of a random effect matrix Z is optional. 

A A 

The only appropriate matrix C for model (3.9) is C ■ [1]; hence where 

a^ is an estimate of the stratum's at-harvest crop acreage proportion for 
year 2. 
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4. DOCUMENTATION OF THE SAS IMPLEMENTATION 


The implementation of the methodology proposed In this paper has been 
developed using the PROC MATRIX feature of SAS (release 79.1) as Implemented 
on the Earth Observations Division Laboratory System (EOOLS) AS/3000 computer. 
It currently resides In accounts JSC1740 and DS40 under the progran Identifier 
GMYP, ''A':. The designation of variable names within the progrwn corresponds ;s 
closely as possible to the names and symbols used In this document. In 
addition, the program has a thorough Internal documentation. 

4.1 REQUIRED INPUTS 

The GMYP/SAS program requires the following matrices as Inputs: PARM, P, X, 

Z, and C. 

The PARM matrix Is a l-by-4 array vector containing the values of operating 
parameters used by the algorithm. Specifically, 

a. PARM(1,1) » the nunber of reweighting Iterations that should be performed 
In arriving at the final estimate of W (usually 2 Is a sufficient value.) 

b. PARM(1,2) » a minerical designator Indicating the transformation that 
should be employed (0, Identity; 1, log; 2, logit). 

c. PARM(1,3) » a specification for the value of cj (usually .001). 

d. PARM(1,4) ■ a specification for the value of C 2 (usually .01). 

The definitions of matrices P (I.e., £), X, Z, and C are the same as those 
given In this document. 

The current version of the program requires that the user Initialize the 
matrices using tlw SAS assignment declarations within the program. Hence, the 
user must be familiar with the procedure for assigning values to matrices In 
PROC MATRIX. The code, however, may easily be modified to accept Input at the 
time of execution. 
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I 


If a model with no random effect is to be fitted, Z should be a column vector 
of zero's (I.e., Z • 0;,.) Also, the program verifies that the specified 
designs satisfy assumptions 1 and j of section 2. If not, processing 
terminates. If the specified design satisfies assumptions 1 and j, process 1m 
will continue even if the resulting model and estimates have no reasonable 



4.2 PROGRAM OUTPUTS 

The GMYP/SAS program provides the user with a mmber of outputs. Including the 
Input data. Intermediate parameter estimates, and final results. All outputs 
are generated using the PRINT option of PROC MATRIX. Hence, all mmerlc quan 
titles are output In matrix format and are labeled by their Interval variable 
names. Table 4-1 shows the variables for which values are output and the 
interpretation of these variables. 



4.3 EXAMPLE RUNS 

A sample input for examples 1, 2, and 3 of section 3 Is shown in figure 4-1. 
Sample outputs for examples 1, 2, and 3 are given In figures 4-2, 4-3, and 4-4. 
The lines of code shown In figure 4-1 should be inserted Into the SAS program 
at the location indicated on the progran listing given In figure 4-5. A flow 
chart showing the functional flow of GYKP/SAS is shown In figure 4-6. 


I 

f; 


4-2 


TABLE 4-1.- VARIABLES OUTPUT IN THE GMYP/SAS PROGRAM 


Variable 

PARM 

P 

X 

Z 

Y 

W 

SIGMAE 

SIGMADE 

GAMMA 

B 

BVAR 

C 

YA 

YAVAR 

A 

ABiAS 

AMSE 

PDIF 


Interpretation 

Row vector of parameters 
n-by-1 Vector ot observed responses 
n-by-t Fixed-effect design matrix 
n-by-s Random-effect design matrix 
Transfornv'd-response vector 
Final iteratively refined estimate of W 
2 

Estimate of a 

e 

2 

Degrees of freedom upon which the estimate of is based 
Estimate of y 
Estimate of ^ ^ 

Estimated V-C matrix for ^ 
c-by-t contrast matrix 
Estimate of C_g = 

A 

Estimated V-C matrix for ^ 

A 

Estimated value of a_ = £ 

A 

Estimated bias of ^ 

Estimated MMS for £ 

Vector of prediction residuals between ^ and its corresponding 
mixed model estimate. (PDIF can be used to calculate measures 
of "goodness of fit" for the model and transformation used.) 
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* DATA FOR EXAMPLEI IN IMPLEMENTATION DOCUMENT* 

PAPM =2 1 *001 •ol; 

PS.279/.154/.149/.07A/.073/.229/.212/.275/.152/.073/.069/.210I 


X=1 

0 

Z=1 

0 

C*1 


0 0/0 
1 0/1 0 
0 0/0 1 
0 1/1 0 
0 0/0 1 


0 

0 

0 

0 

0 


0/0 0 1 
1/0 1 0 
0/0 1 0 
0/0 I 0 
0/0 0 1 


0/1 

1/1 

0/0 

0/0 


0 0 0/0 
0 0 1/0 
0 1 0/0 
0 1 0/0 


100/01 

101/00 


0/0 

0/0 


0/1 -1 00/001 i ; 


0 0 / 
1 1 * 


0 

0 


1 / 

It 


a. Sample input, example 1. 


* DATA FOR EXAMPLE2 IN IMPLEMENTATION DOCUMENT* 
PA»M=2 1 .001 .01? 


Ps. 

279/. 

15^/. 1^9/. 074/ 

• 

073/ 

.229/ 

.212/. 275/. 152/. 

073/ 

.069/. 210* 

X=1 

0 

0 

0 1 0/01 

0 

0 

1 

0/0 

010 1 0/1 

0 

0 0 

-1 

0/ 

0 

1 

0 

0 -1 0/0 1 

0 

0 

-1 

0/0 

010-1 0/1 

0 

9 1 

1 

1/ 

0 

1 

0 

1 1 1/10 

0 

1 

-1 

-1/0 

101-1 -1/0 

0 

1 1 

-1 

-1* 

Z=1 

0 

0 

0/0 1 0 0/0 

1 

0 

0/0 

0 1 

0/0 010/00 

0 

1/ 



0 

0 

0 

1/1 0 0 0/0 

1 

0 

0/0 

0 1 

0/0 010/00 

0 

1* 



Csl 

0 

0 

000/010 

0 

1 

0/0 

0 1 

0 -1 0* 






b. Sample input, example 2. 


* DATA FOR EXAPMLE3 IN IMPLEMENTATION DOCUMENT* 
PARM=2 2 .001 .01* 

P=.15A/.073/.229* 

X=l/1/1 t 
Z=0/0/0l 
C=l* 


c. Sample input, example 3. 


Figure 4-1.- Sample input for examples 1, 2, and 3 of section 3. 
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PARM 

COLl 

C0L2 

COL3 

C0L4 

ROwl 

2 

1 

0.001 

0.01 

P 

COLl 




RUWl 

0.279 




ROw2 

0.1S4 




ROw3 

0.149 




ROw4 

0.074 




ROWS 

0.073 




R0W6 

0.229 




ROw7 

0.|l2 




ROWS 

0.275 




ROw9 

0.152 




ROmIO 

0.073 




ROwU 

0.069 




R0wl2 

0.21 





X 

COLl 

C0L2 

COL3 

COL4 

ROwl 

1 

0 

0 

0 

ROw2 

0 

1 

0 

0 

R0W3 

0 

0 

1 

0 

ROw<» 

1 

0 

0 

0 

ROWS 

RUw6 

0 

0 

i 

0 

0 

0 

0 

ROW? 

0 

0 

1 

0 

ROWS 

1 

0 

0 

1 

ROw9 

0 

1 

0 

1 

ROwlU 

1 

0 

0 

i 

ROW 11 

0 

1 

0 


R0wl2 

0 

0 

1 

i 


z 

COLl 

COL 2 

COL 3 

COL4 

ROwl 

1 

0 

0 

0 

ROW2 

0 

1 

0 

0 

R0W3 

0 

1 

0 

0 

ROWh 

0 

0 

1 

0 

ROWS 

0 

0 

I 

0 

ROWb 

0 

0 

0 

I 

ROW? 

0 

0 

0 

1 

0 

0 

ROwB 

1 

0 

0 

R0W9 

0 

1 

0 

ROW 10 

0 

0 

1 

0 

ROWll 

0 

0 

1 

0 

ROW12 

0 

0 

0 

1 


RO«l 

RUw2 

ROn3 

RO«<» 

ROWS 

ROWb 

ROW? 

R0W4 

R0W9 
ROwlO 
ROwIi 
ROW 12 


COLl 

•1.2765A 

-1.8708 

•l.90?81 

■2.60^69 

-2.6173 

•l.A7«»03 

’i2?ll73 

Figure 4-2.- Sample output 'or example 1. 
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COLl 

COLT 

S§ti 

SOwl 

0.368701 

0 


U 

0 

R0m2 

0 

0.184369 


0 

0 

RUm3 

0 

0 


0 

0 

ROw<» 

0 

0 


0 

0 

R04S 

0 

0 


0 

0 

ROm6 

0 

0 


0 

0 

R0«7 

0 

0 


0.27259 

0 

ROWS 

0 

0 


0 

0.377576 

R0m9 

0 

0 


0 

0 

R0«1U 

0 

0 


0 

0 

ROW 11 

0 

u 


0 

0 

ROW12 

0 

0 


0 

0 


IS SYSTEM 


C0L3 
COL 9 

C0L4 

COLIO 

C0L5 

COLll 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.17213 

0 

0 

0 

0 

0 

0 

0.0802295 

0 

0 

0 

0 

0 

0 

0.0772788 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

oo 

oo 

0 

0 

0 

0 

0 

0 

0 

0 

u 

0.179649 

0 

0 

0 

0 

0.07B431§ 

0 

0 

0 

0 

0 

0 

0 

0.0755519 

0 

0 

0 

0 

0 

0 


5I6MAE 

COLl 

ROwl 

0.000033604 

SIGHAOr 

COLl 

ROWl 

5 

GAMHA 

COLl 

ROWl 

10992.7 




0 

0.abS<»23 





Figure 4-2,- Continued. I 

•I 
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BVAR 

COLl 

COL2 

ROWl 
ROW 2 
ROw3 
ROw4 

0.0925477 

0.0922792 

0.092255H 

-.000022708 

0.0922792 

0.0924419 

0.0924181 

-.000022697 

C 

CULl 

COL2 

row! 

RUw2 

ROn3 

RUW4 

ROWS 

1 

0 

0 

1 

0 

0 

1 

0 

-1 

0 

YA 

COLl 


ROwl 

R0W2 

ROW3 

ROW4 

ROWS 

-1.78651 

0.0348533 

-1.90061 


YAVAR 

COLl 

COL2 

ROwl 

R0w2 

RUw3 

ROw^» 

RUw5 

0.0925477 

0.0922792 

0.0922554 

0.000263564 

0.0922327 

0.0922792 

0.0924419 

0.0924181 

-.000162763 

0.0923954 

A 

COLl 


ROWl 

R0W2 

ROw3 

RUW4 

ROWS 

0.167543 

0.161805 

0.15265 

1.03547 

0.149477 


ABIAS 

COLl 


ROwl 

ROw2 

ROW3 

ROWi. 

ROWS 

0.00775288 

0.00747876 

0.00706043 


AMSE 

COLl 

COL2 


ROnI 

R0k2 

ROM3 

ROW4 

ROi^S 


C0L3 

0.09225134 -.OOOOi 
0.0924181 *.0U00< 


COL3 

0 

0 

1 

0 

1 


C0L4 

illSf 


C0L4 

0 

0 

0 

0 

1 


COL3 


COL4 


0.0922554 0.000268564 
0.0924181 -.000162763 


0.0924749 -.000162726 


COL5 

0.0922327 
0.0923954 
0.0924749 
«. 00016272b 
0.0925046 


0.00259789 

,0000465919 

0.00230987 


0.00250162 

0.00242Q19 

0.0022326t 


C0L3 

0.00235948 
0.00228268 -( 
o.ooSif 


5556 -( 


COL4 

000465919 
‘02727 
02572 


•O.OOOQ2727 -0.00002572 0.000462466 
0.00223468 0.00211007 -.000025187 


C0L5 

0.00230987 
0.0Q223466 
0. 00211007 
-.000025187 
0.00206687 


Figure 4-2.- Continued. 


4-7 


ORIGINAL PAGE IS 

OF Fv)OR QUAUT/ 


statistical analysis system 


PDIF 

POWl 

P0m2 

ROn3 

ROH4 

RUii«5 

ROWS 

R0W7 

ROwB 

R0W9 

ROW 10 

Rowii 

ROW12 


COLl 

-oJoiUols 

-0.0036504S 

-0.093SA33 

-0.0888045 

0.0671955 

0.0593495 

0.110939 

-0,00644(096 

-0.0910605 

-0.089441 

0.0605228 


Figure 4-2.- Concluded. 
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MAMH 


COLl 


COLi! COl3 CUL<* 


RUWl 




O.OUl O.Ol 


p 

COLl 


ROwl 

0.279 


R0W2 

0. ISA 


ROW3 

0. 1A9 


ROw4 

0.07<» 


ROWS 

0.073 


ROWS 

0.229 


ROW7 

0.212 


ROWS 

0.27S 


R0w9 

0. 1S2 


ROwl 0 

0.073 


ROWl 1 

0 . Onv 


ROW12 

0.21 


X 

COLl 

COL? 


COl3 cula cols cols 


ROwi 

WUW2 

RU«3 

HUW4 

f^O «' b 
RUW6 
ROw7 
ROWS 
RUWR 
ROWK! 
KO«l ) 
ROW 12 


1 

0 

0 

1 

u 

0 

0 

1 

0 

I 

0 

0 


0 

1 

0 

u 

1 

I 

0 

0 

1 

u 

1 

0 


0 

0 

I 

J 

u 

0 

1 

0 

u 

0 

0 

1 


0 

0 

0 

u 

0 

0 

0 

I 

1 

1 

1 

1 


1 

1 

1 

-1 

-1 

“1 

-I 


-1 

-1 

-1 


1 

1 

-1 

-1 

-1 


z 


COLl COL2 


COL3 COLa 


ROWl 

ROW2 

ROW3 

ROW A 

ROWS 

R0W6 

HOW7 

ROWS 

ROWO 

ROW 10 

ROWII 


1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 


u 

I 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 


0 

u 

0 

1 

1 

0 

0 

0 

0 

1 

1 

0 


0 

0 

0 

0 

0 

I 

0 

0 

0 

0 

1 


y 

COLl 

ROWI 

-1.2765A 

R0w2 

-1.8708 

ROW3 

-1.90361 

ROWA 

''2.60 369 

ROWS 

'2.6173 

ROWS 

.A7A03 

ROW? 

.SSI 17 

ROWS 

. .29098 

R0W9 

-1.88387 

ROW 10 

-2.6173 

row! 1 

-2.67365 

ROW 12 

-1 .S606S 


Figure 4-3.- Sample output for example 2. 
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w 

con 

C0L2 

COL3 

CULW 

CuLb 

.C0L6 


COL? 

COLb 

COLV 

COLlU 

COL 11 

COL 12 

RO*(l 

0.387995 

0 

0 

u 

0 

0 

0 

U 

0. 

u 

0 

0 

R0W2 

0 

0. lB40b7 

0 

u 

0 

0 

0 

0 

0 

0 

0 

0 

ROW3 

0 

0 

0.172151 

u 

0 

0 

0 

0 

u 

0 

0 

0 

R0M4 

0 

0 

0 

0.0803?SA 

0 

0 


0 

0 

u 

0 

0 

0 

ROWS 

0 ’ 

0 

0 

0 

0.0774223 

0 

0 

0 

u 

0 

0 

0 

RUW6 

0 

0 

0 

0 

0 

0.29379b 


0 

0 

0 

0 

0 

0 

ROW? 

0 

0 

0 

u 

0 

0 


0.273154 

0 

0 

0 

0 

0 

ROWe) 

0 

0 

0 

0 

0 

0 

0 

0.3?a2d2 

0 

0 

0 

0 

HOW 9 

0 

0 

0 

0 

0 

0 


0 

0 

0. 1h0122 

0 

0 

0 

ROWIO 

0 

0 

u 

u 

0 

0 

0 

0 

u 

0.0/b2«<»tt 

0 

0 

ROW 11 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0. 0754138 

1) 

R0W12 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.2b4bl9 


SI6MAE 

COLl 


ROWl 

.0000408873 


S16MADF 

COLl 


ROWl 

4 


gamma 

COLl 


ROWl 

15874.5 


B 

COLl 


ROWl 
ROW2 
ROW3 
ROW 4 
ROWS 
ROW6 

“1.78686 
-1.82155 
-1.87823 
-0. 0213883 
0.23617 
0.0030828 

Figure 4-3.- Continued. 
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HVAP 

cou 

COL^ 

COL3 

cola CQL5 

ROWl 

o.iea'-ort 

0. I6^1H2 

0.1b?lA9 

-.00u0^obAb -.OOOOAJlbT 

HO*<2 

0. 

0. 1!>2 Jtt 

0. le>23'*7 

-.U00U2bM3/ .00002l?^<'5 

«UW3 

0. lt f’1.9 

0.1()2-H-t7 

0. 7r 

-.U0t)0<.0:>bV .00001b‘.772 

MOWA 

-.oooo?t«fl-.t) 

-.OOOtVbdJ/ 

-.OOOOAObb'i 

.0000 / Jc. 227 .OOOU02h;>AB 

MOWS 

0000*. J1S7 

• 000U212‘.2b 

.OOOOln<*7 Id 

.U00u02Mb**o 0.i(s2.10b 

HUM6 

-6.39{Jbt-06 

-6.A068L-0b 

O.OOUC Jt>26b 

-.U00010H7b -.000030988 


C 

COLl 

C0L2 

C0L3 

COLA 

COLS 

R0*(1 

1 

0 

0 

0 

0 

H0W2 

0 

1 

U 

U 

1 

H0ip<3 

0 

0 

1 

0 

-1 

YA 

COLl 





ROWl 

-1 . 





R0W2 

“1.88838 





RUW3 

“2. llAA 





YAVAH 

COLl 

C0L2 

C0L3 



ROWI 

0 « 1 b f 8 0 8 

0 . 16 2138 

0.162192 



H0W2 

0. lb213H 

0.32A/2M 

.00003633'^ri 



R0W3 

0.162192 

.0000363368 

0 . 32*» 7 ‘**» 




A 

COLl 

ROWl 

0.167485 

H0X2 

0.2048/ 

R0M3 

0.120705 


ABIAS 

COLl 

ROWl 

0.0136068 

RUW2 

0.033263b 

ROW 3 

0.0195992 


AMSE 

COLl 

C0L2 

COl3 

ROWl 

0.00455655 

0.0055634 

0.00327893 

8.9857t-07 

ROW2 

0.0055634 

0.0136295 

ROW3 

0.00327993 

6.9857t-07 

0.00A73145 


PDIF 

COLl 

ROWl 

0.066699 

R0W2 

-0.0508705 

H0W3 

“0.0445HU2 

R0W4 

“0*0582538 

ROWS 

-0.05A7453 

RUW6 

0.101255 

ROW? 

0.0912947 

ROwe 

0.j6b74b3 

ROW9 

-0.0<.91543 

ROWIO 

-0.0560567 

row! 1 

-0.0556572 

ROW 12 

0.092212b 


Figure 4-3.- Concluded. 


COLb 

fe.a^JrtSE-Ob 
b.'.OhoF-Ob 
. OOOO 

.OOOUIOh 7b 
. UUUO J09H6 
.00UUU« 107 


COL6 

0 

0 

0 
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b T A T 1 

1 b T I C A L 

A U A L 

PAKM 

COLl 

C0L2 

KUWl 

2 

2 

P 

COLl 


POWl 

0.1b'. 


POW2 

0.07 J 


POW3 

0.22'y 


X 

COLl 


KUnI 

POW2 

KOW3 

1 

1 


2 

COLl 


MOWl 

0 


POW2 

0 


ROn3 

0 


y 

COLl 


MOW! 

-0.6bl7w3 


«OW2 

-1.2707b 


ROW3 

-0.60b9H3 


w 

COLl 

C«JL2 

ROWl 

0.12002 

0 

MOW? 

0 

0.12002 

MOW 3 

0 

0 

SIGMAE 

COLl 


ROW] 

0.0135231 


SIGMADF 

COLl 


ROW! 

2 


gamma 

COLl 


ROWl 

0 


B 

COLl 


MOWl 

-0.909838 


BVAR 

COLl 


ROW] 

0.0375578 



lb S f b T t M 
COLJ 

0.001 O.Ol 


C0l3 


0 

0 

0.1200<f 


Figure 4-4.- Sample output for example 3. 
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c 

COL) 

ROW) 

1 

YA 

COLl 

ROW) 

-0.909P3tt 

YAVAR 

COLl 

ROW) 

0.0375578 

A 

COLl 

HOW) 

0.139A73 

A8IAS 

COLl 

ROW) 

0.00650056 

AMSE 

COLl 

ROW) 

0.00216A05 

PDIF 

COLl 

ROW) 

R0W2 

ROW3 

0.01A5272 

“0.066A728 

0.0895272 


wage iS 
quality 


Figure 4-4,- Concluded. 
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lOPAlNG VERSION - A*PlW rMcrn ObAlCiOit 

•HOGHAMMEO BY! T.BAKtw / LOCKMEEO-EMSCO <Za5I69)I 




• SAS PPOGPAM TO 


Y(PI 


»»fl • 2*R 


El 


•I 





proportion ESTIMATES! 


OR OF F 


lA F 

ISn PftTRi5^?OM\'THC_§ANOOM EF 


EfrecTst 


♦I 


# 

# 




•I 


model tfTIiTliiooM vector -ith, 
- ?4io5Ti^ONAJj^T5J|.(l4j.! 


2) 

3) 

A) 

5) 


6 ) 

7) 

B) 


MEAN PI AND VARIANCE! 
Ih' mean ZERO AND DIAGONAL 


R Ts A RANOOM VEClun ■ ! ' t TY . - 

idc?^ifSivsIolr§^^ 

' r MODEL PITH NO RANDOM EFPeCT IR) IS TO 6E PlfUDM 
IhE NUMBER OF *08SERV at IONS IS AT L^ST RANK ( X : 21 • 1 1 


C MATRIX! 
MATRIX! 


REQUIRED INPUTS! PARMf P 
OUTPUTS: 


X< Z* C! 


PAHMt Pt X* 2* 7» R* SI6MAE. 
Ct YA< YAVAR* At AblASt AMSC* 


^^E7fT' 


gamma t Bt BVAR! 


internal VARIAB 

PARM b a 

PAHM ( 


•‘hS-'^VElforSf^StEPllHM PARAMETERS! 

. 1 ) !n:mbEP of RE veTghT I NG I Tt« atTonS 

, 5) bTHANSFORMATION.TyPE (O-IUEnUTy 1 1 


PARM( 

PAPM I 
PARM( 

■ THE 

B TmE 

B TmE - 

A matrix rmo 


(USUALLY^?) 
y06t4-LO$l 


I,! 


i^iiToJtRANCE^LlMlT'^ofe .Q9 

.AliTOLERtNcI lImIt for RElGnt CalCULAtIgnI (USUALLY .51 
COLUMN VECTOR OF HESPOnSESI 
FIXED EFFECTS DESIGN MATRIX! 

R ANUOM effects. 0E|l6N,M*TR I * 


)l) ! 

U ! 


^Ovs‘^^ON§15T"6F*FixED EFFECTS DESIGN SETTINGS! 
FOR_RMi^M TRANSFORmIo-RESPOnIe IYA) AND! 


RESPONSE (A) 

^ VECTOR! 

- A PORK matrix hHQS 
matrix FOR THE M'*''-* 

IGmaE ■ THE ESTIm 

issr. i-Ii|l« 2 feS* 2 ilSS 6 = 

“?ir5Ic?5h;.ri.554rr|jT5STj-jT|f. 


iKi 5 s»ib_? 85 iiyi|y 5 ! 


THE LAST UPDATED BEIOhTING! 


lOR VARIANCE 
^F FREEDOM! 


0 TRAN 


VAM(SORTIWI«E) ! 
UMMED'ERROR! 


YAVAR » THE ESTIMATtO VAKia 

A i**rMpyJcT0H*0F INVERSEP-TRANSFORmE! 

s the PREOlCTtD RESPONSES FOR THE 0 
ABIAS » A VtcTOR OF APPROXIMATE 


MATRIX FOR Bl^^ 

(Y-VALUES) FOR! 

X^FOR The PREDICTED! 


* matrix for 


PDIF B P - PHAT B 


NR 

NRC 


number of ROMS (OBSERVATIONS) 

... ■ NUMBER “■ '• 

NCX s NUMBER 

Ncz B number . .. 

X 2 B THE MNCATtNA. 
rWx 2 b Rank of (Xi 2 ) ••(A! 2 | ! 
PHAT ■ A w^wK VECTOR BHOSE r I 


C! 

ANCE'maTrIx for the PREDICTED! 
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Figure 4-5.- Program listing for program GMYP/SAS. 
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Figure 4-5.- Continued. 
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Figure 4-5.- Continued. 
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Figure 4-5.- Concluded. 
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Figure 4-6.- Functional flow chart for program GMYP/SAS. 
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